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| We report numerical and analytic results for the spatial survival probability for fluctuating one- 

. dimensional interfaces with Edwards- Wilkinson or Kardar-Parisi-Zhang dynamics in the steady 

state. Our numerical results are obtained from analysis of steady-state profiles generated by in- 
tegrating a spatially discretized form of the Edwards- Wilkinson equation to long times. We show 
that the survival probability exhibits scaling behavior in its dependence on the system size and 
«4 ■ the "sampling interval" used in the measurement for both "steady-state" and "finite" initial con- 

ditions. Analytic results for the scaling functions are obtained from a path-integral treatment of a 
formulation of the problem in terms of one-dimensional Brownian motion. A "deterministic approx- 
imation" is used to obtain closed-form expressions for survival probabilities from the formally exact 
analytic treatment. The resulting approximate analytic results provide a fairly good description of 
C$ , the numerical data. 
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S ■ I. INTRODUCTION 

Temporal first-passage properties, expressed in terms of persistence and survival probabilities 1], have recently 
found many applications QUO, 01 m * ne study of the dynamics of fluctuating interfaces. Experimental realizations of 
one-dimensional (Id) fluctuating interfaces are provided by monatomic steps on vicinal surfaces. Recent experimental 
studies 0, 0, H, El have shown that temporal persistence and survival probabilities for fluctuating steps can be 
measured, and that these probabilities provide a convenient way of characterizing their dynamics. The temporal 
persistence probability is defined in this context as the probability that the height at a particular point of the 
interface does not return to its initial value over a certain period of time. For many simple models of interface 
dynamics, this probability exhibits a power-law decay in time [2j.|3|.Ll|. This power-law behavior has been confirmed 
in experiments [E 0, IS 13 In contrast, a closely related quantity, the temporal survival probability that measures 
the probability that the height does not return to its average value over a certain period of time, is found, both 
lO [ theoretically and experimentally H,|I1 to decay exponentially at long times. 

In studies of fluctuating interfaces, it is natural to consider spatial analogs of these temporal first-passage quantities, 
| namely the spatial persistence and survival probabilities. These probabilities have been studied analytically |l(J 

S and numerically [jjj for several models of interfacial dynamics. For (l-t-l)-dimensional interfaces, the stochastic 
variable of interest is the "height" h(x, t) that represents the position of the interface at point x and time t. In this 
. paper, we consider the interface profile h(x,to) where the time to is in the long-time, steady-state regime, and for 
notational convenience, suppress the time argument to of h from now on. To define spatial persistence probabilities, let 
p(xo, xq + x) be the probability that the height h(x) does not return to its "initial" value h(xo) at the point xo within 
O ■ a distance x measured from xq along the interface. The "steady-state" spatial persistence probability P(x) is then 
defined as the average of p(xq, xq + x) over all possible choices of the initial point xq in a steady-state configuration 
of the interface. A seco nd p ersistence probability Ppjc(x), the so-called spatial persistence probability for "finite 
initial conditions" (FIC) is defined as the average of p(xo, xq + x) over initial points xq at which both the height 
h(xo) and its spatial derivative h'(xo) are finite. It was shown in Refs. [ToL ITlj that for several models of fluctuating 
interfaces, both P(x) and Pfic(x) exhibit power-law decay for large x, but the exponents that describe this power- 
law behavior may be different in the two cases. Experimental measurements of the steady-state spatial persistence 
probability for interfaces (combustion fronts in paper) believed to be described by the Kardar-Parisi-Zhang (KPZ) 
equation [l^ have been reported recently jl3j |. The behavior of Pfic( x ) wa s not investigated in this work. 

In analogy with the temporal case, the spatial survival probabilities are defined in terms of the probability p' (xq , xq + 
x) that the interface height between points xq and Xq + x does not cross its average value h (rather than the initial 
value h(xo))- The steady-state and FIC spatial survival probabilities, Q{x) and Qfic(x), respectively, are then 
obtained by averaging p'(xq, xq + x) over xq in the two different ways mentioned above: in the first case, the average 
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is done over all points Xq, while in the second case, the average is performed over only the points at which the height 
and its spatial derivative are finite. Numerical results for the x-dependence of these two spatial survival probabilities 
for Id interfaces with KPZ and Edwards- Wilkinson (EW) dynamics were reported in Ref.^3|. It was found 
there that the spatial dependence of Q(x) is neither power-law, nor exponential, while Qfic{%) exhibits a power-law 
decay similar to that of the FIC persistence probability Pfic(x). While the power-law behavior of Qfic{x) was 
expected the rr-dependence of Q(x) found in the numerical work of Ref.^]| was not understood theoretically. 

In this paper, we present the results of a detailed numerical and analytic study of the spatial survival probabilities 
for Id interfaces with EW dynamics in the steady state. The primary motivation for this study is to develop a 
theoretical understanding of the numerical results reported in Ref. 01 ■ On a more general level, studies of temporal 
and spatial first-passage properties of fluctuating interfaces are believed to be important in understanding the role 
of thermal fluctuations of the edges of components in the stability of nano-scale devices. The Id EW equation is 
believed |l5j to describe thermal fluctuations of steps on a vicinal surface under experimental conditions for which 
the dominant source of fluctuations is the attachment/detachment of atoms to/from the steps. Since the statistics 
of height fluctuations in the steady state of the Id KPZ equation is the same as that for the Id EW equation, our 
results also apply to experimental realizations of Id KPZ interfaces. 

The Id EW equation has the form[l4j 

^M = r ^M + ri (x,t), (i) 

where T is a kinetic parameter and rj(x, t) is a Gaussian random noise with (r)(x, t)) = 0, (r)(x, t)r)(x' , t')) — 2D' 5{x — 
x')6(t — t'), D' being a measure of the strength of the noise. The parameters T and D' should satisfy the fluctuation- 
dissipation relation if this equation is supposed to describe equilibrium fluctuations (e.g. in the case of steps on a 
vicinal surface), but they are independent parameters in a general non-equilibrium situation. We consider a finite 
system of size L, so that < x < L, and use periodic boundary conditions. It is easy to see that the spatial average 
h(t) = l/L J Q h{x,t)dx executes a simple random walk in time. In our calculations, we subtract the spatial average 
from the variables h[x,t), so that from now on, it is implied that h(x,t) represents the height at point x measured 
from the instantaneous spatial average. Thus, h(t) is equal to zero by definition. As we shall sec later, this condition 
plays an important role in the analytic calculation. 

The spatial survival probabilities studied here are defined as follows. Let Q{x, L\ho) denote the probability that 
the steady state spatial profile of the interface starting at h$ at x — (and ending at h$ at x = L, i.e., with periodic 
boundary condition) does not cross upto a distance x where < x < L. Then, the steady-state spatial survival 
probability is given by 

/oo 
Q(x,L\h Q )P st (h ,L)dh , (2) 
-OO 

where P st (ho,L) is the steady-state height distribution. The FIC survival probability is defined as 

n . T , f™ w Q(x,L\h )P st (ho,L)dh 

Qfic{x,L,w) = w — — — , (3) 

J- w Pst{ho,L)dh 

where w <C W{L), the steady-state width of the interface. In our numerical work, steady-state profiles for systems of 
different L are generated by integrating a spatially discretized form of the EW equation to long times, and these profiles 
are used to calculate the survival probabilities. The spatial discretization implies that there is a finite sampling interval 
8x that represents the spacing between successive points at which the height variable is sampled in the calculation of 
the survival probabilities. Clearly, the value of Sx must be an integral multiple of the spacing of the spatial grid at 
which the height variable is defined. This sampling interval Sx is analogous to the "sampling time" that represents 
the interval between two successive measurements of a stochastic process in studies of temporal persistence. The fact 
that a finite value of the "sampling time" may modify the persistence properties of a stationary stochastic process 
was first pointed out in Ref. fl6| . In the context of fluctuating interfaces, it is known from numerical 0,0] and 
experimental [9j studies that the temporal survival probability in interfaces exhibits a non-trivial dependence on the 
sampling time. In our numerical study, we find that the spatial survival probabilities also depend on the value of 
the sampling interval Sx. The dependence of the survival probability Q on x, L and Sx is found to be described by 
a scaling function of x/L and Sx/L: Q{x,L,5x) — fd{x/L,8x/L). This is similar to the scaling behavior found in 
Ref.[lll| for spatial persistence probabilities. As shown there, the dependence of the FIC survival probability Qfic 
on x, L, Sx and w is also described by a scaling function of x/L, Sx/L, and w/L a where a = 0.5 is the exponent for 
the dependence of the steady-state width W(L) on L [W(L) cx L a }. 
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In our analytic study, we consider the spatial survival probabilities when the height variable is sampled continuously, 
i.e. the limit Sx — > 0, and calculate the scaling function f{x/L) = fd(x/L, 0). This calculation is based on a mapping 
of the spatial statistics of a steady-state EW interface to the temporal statistics of Id Brownian motion. The 
requirement that the average height h of the interface must vanish translates in this mapping to the constraint that 
the total area under the curve that represents the Brownian process in the distance - time plane must be zero. This 
"zero-area" constraint plays a very important role in the analytic calculation - the form of the scaling function f(u) 
depends crucially on whether this constraint is imposed in the analytic treatment. Without this constraint, we can 
determine f(u) exactly, but the form of the scaling function obtained this way differs drastically from the numerical 
result. In particular, the scaling function does not go to zero as u — > 1, whereas the numerically obtained scaling 
function decreases rather fast to as u approaches 1. When we take into account the zero-area global constraint, 
determining the scaling function f(u) analytically becomes much more nontrivial. We are able to set up an exact path 
integral technique that allows us, in principle, to compute this function exactly in terms of some complicated integrals. 
However, we cannot get an exact closed form expression for f(u) to compare with the simulation data. We then make 
a simple "deterministic" approximation that allows us to obtain a a closed form expression for f(u) which we then 
compare with the numerically obtained scaling function. The agreement is fairly good, given the drastic nature of 
the deterministic approximation. Our approximate analytic result for the FIC survival probability also shows similar 
agreement with the numerically obtained results. 

The rest of this paper is organized as follows. In Scc.[H] we describe our numerical results for the spatial survival 
probabilities. The analytic calculations with and without the "zero-area" constraint are described in detail in Sec. lIIII 
In this section, we also present a comparison of the analytic results for the survival probabilities with the numerical 
results presented in Sec. [n] Sec. IIVI contains a summary of the main results and a few concluding remarks. Some 
details of the steady-state properties of finite Id EW interfaces with periodic boundary conditions are presented in 
the Appendix. 



II. NUMERICAL RESULTS 



In the numerical work, we consider a spatially discretized dimensionless form of the Id EW equation defined in 
Eq. l[T]l. The height variable is defined on a Id lattice of unit spacing with periodic boundary conditions. Let hi be 
the height at lattice site % with i — 1, 2, . . . , L. The time-dependence of the height variables is given by 

^1 = [h i+1 (t)-2h i (t)+h i . 1 (t)]+r ]i (t), (4) 
at 

where the r/i(t)'s represent uncorrelated Gaussian noise with (r)i(t)) = and (%(t)?7j(t')) = 2<5y<5(i — t'). The Eq. Q 
is thus a discretized version of the continuum EW equation Q| w ith the choice T = 1 and D' = 1. These equations 
are integrated forward in time using the simple Euler method |17| . Thus, we write Eq. Q as 

hi(t + St) - hi(t) = 5t[h l+1 {t) - 2hi(t) + hi-!(t)] +Vsin(t), (5) 

where each rj(i) is an independent Gaussian random number of zero mean and variance equal to 2. We use a value 
of St small enough (i.e. St = 0.01) to avoid any numerical instability. Steady-state interface profiles are generated 
by carrying out the integration from flat initial states (hi — for all i) to sufficiently long times (much longer than 
the time at which the width of the interface saturates). As mentioned in Sec[IJ we always subtract the instantaneous 
spatial average of the height variables from the individual variables {hi}, so that the condition XaLi = is always 
satisfied. 

The steady-state spatial survival probability Q(x) is measured as the probability that the interface height variable 
does not cross zero as one moves along the interface from an initial point Xq to the point xq + x (since the height 
variables are defined on a lattice of unit spacing, both xq and x are integers between and L.) This probability 
is averaged over all initial points xq in a steady-state configuration and also over many (typically 10 4 ) independent 
realizations of the stochastic evolution that generates the steady-state profiles. The minimum value of the sampling 
interval Sx used in the measurement of Q(x) is obviously the lattice spacing which is equal to unity. However, it is 
also possible to use a larger sampling interval, equal to a positive integer m, in the measurement of Q(x) - this is 
done by considering only the heights at the lattice sites i = km, k — 1,2,... while checking whether the height crosses 
zero between the points xq and xq + x. The measured survival probability exhibits a weak dependence on the value 
of Sx. The FIC survival probability is measured in a similar way, except that the initial points xq are chosen to be 
only those at which the height lies between — w and +w, with w <C W(L), the steady-state width of the interface. 

Typical results for the survival probability Q(x) are shown in Figs. 11141 In Fig. 2] we show plots of Q{x) as a 
function of x/L for L — 200, 400 and 800. The values of Q{x) shown in these plots were obtained using Sx = 1 (unless 
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FIG. 1: (Color online) Plots of the spatial survival probability Q(x,L) as a function of x/L (L is the sample size) for L = 200 
((black) dotted line), L = 400 ((blue) dash-dotted line), L = 800 ((magenta) dashed line), and L = 10 4 ((red) full line). The 
same sampling interval, Sx = 1, is used in all cases. 



mentioned otherwise, all the results shown here were obtained using this "default" value of Sx). It is clear from the 
plots that Q(x) decreases from 1 to a value close to as x/L is increased from to about 0.6. The numerical results 
show a weak dependence on the value of L. As we shall see shortly, this dependence arises from the use of the same 
5x in all the measurements for the different values of L. 

The numerical calculations can not be extended to much larger values of L because the time required to reach 
the steady state from a fiat initial state increases with L as L z with z = 2. However, we have found a different 
way of generating steady-state interface profiles for much larger values of L. It is easy to show that in the steady 
state of the model defined in Eq. the height difference variables Sj = h i+ i — hi, i = 1,2, ... ,L (with h^+i = hi 
due to periodic boundary conditions) are independent Gaussian random variables with zero mean and variance equal 
to unity, apart from the obvious constraint, X)i=i s « = 0, arising from periodic boundary conditions. Therefore, a 
realization of the steady-state interface profile for a system of size L may be obtained by numerically generating L 
Gaussian random variables with the statistics mentioned above, identifying these random variables with the s,'s, and 
then calculating the heights /i, (with their spatial average subtracted off) in terms of these Si's. We have calculated 
the spatial survival probabilities for steady-state EW interfaces with L = 10 4 generated this way, averaging over 4000 
independent realizations. The results for Q(x) obtained from this calculation are also shown in Fig^ It is clear 
that these results for L = 10 4 are consistent with the trend shown by the other results obtained from steady-state 
interfaces generated by numerical integration. 

In FigO we have shown plots of Q(x) for L = 800, obtained from two different calculations, one with Sx = 1 and 
the other with Sx = 4. The two curves are clearly different, indicating that the measured Q(x) depends weakly on 
the value of Sx used in the measurement. We have found, in analogy with the results reported in Ref.|ll| for spatial 
persistence probabilities, that the dependence of Q on x, L, and Sx satisfies the scaling equation 

Q(x, L, Sx) = fd(x/L, Sx/L), (6) 

where the subscript "d" of the scaling function is meant to indicate that here we are considering survival probabilities 
measured using discrete sampling with a finite sampling interval Sx. This scaling equation implies that plots of Q(x) 
vs. x/L for samples with different L would all collapse to the same curve if the survival probabilities for different 
L are measured using different values of Sx, such that Sx/L is held constant. This scaling behavior is illustrated in 
FigEl I n this Figure, we have shown plots of Q vs. x/L for L = 200 measured with Sx = 1, L = 400 measured with 
Sx = 2, and L — 800 measured with Sx = 4 (so that Sx/L = 1/200 in all 3 cases). The 3 sets of data are found to 
collapse to the same curve, thereby establishing the validity of the scaling form of Eq. © . The L-dependence of the 
results shown in Fig^may, therefore, be thought of as representing the dependence of the scaling function fd on the 
value of its second argument, Sx/L. 
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FIG. 2: (Color online) Dependence of the spatial survival probability on the sampling interval Sx used in its measurement. 
Plots of Q(x,L) versus x/L (L is the sample size) are shown for L = 800, and two values of 5x: 8x — 1 ((black) full line) and 
Sx — 4 ((red) dashed line). 
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FIG. 3: (Color online) The scaling behavior (Eq. JSJ) of the dependence of the spatial survival probability on the sample size 
L and the sampling interval 8x. Plots of Q(x,L) versus x/L are shown for three different sets of values of L and 8x with 5x/L 
held constant: L = 200 and Sx — 1 ((blue) crosses); L = 400 and Sx = 2 ((black) line); L — 800 and Sx — 4 ((red) circles). 

Numerical results for the FIC survival probability Qfic were reported in Ref.^lJ where it was shown that it 
exhibits the following scaling behavior: 

Qfic{x, L, Sx, w) = f F ic{x/L, Sx/L, w/W(L)). (7) 

For the sake of completeness, we have shown in FigQJour numerical results for Qfic obtained for w/W(L) — 0.02. 
Two sets of data are shown, one for L = 800, obtained from steady-state interfaces generated by numerical integration, 
and the other for L = 10 4 , obtained from interfaces generated using Gaussian random variables, as outlined above. 
In both cases, the initial decay of Qfic can be well-represented by a power law, Qfic(%) oc x~ 6fic with Qfic = 0-5. 
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FIG. 4: (Color online) Double-log plots of the FIC spatial survival probability Qfic(x, L) versus x/L (L is the sample size) for 
L = 800 (crosses) and L = 10 4 (circles). Fits of the initial decay of Qfic to a power-law with exponent 1/2 (Qfic(%) oc x~ 0,5 ) 
are shown by the (red) solid lines. 

As in the case of the steady-state survival probability, the dependence of Q pic on the value of L may be thought of 
as representing the dependence of the scaling function /fic of Eq. J7J on the argument 5x/L. 

III. ANALYTIC CALCULATIONS 

In this section, we describe in detail our analytic calculation of the spatial survival probabilities in the stationary 
state of the Id EW equation QJ. We consider periodic boundary condition, h(0,t) = h(L,t). The height can then 
be decomposed into a Fourier series, h(x, t) = ^ fc h(k, t)e tkx where k = 2nm/L with m — 0, ±1, ±2 .... Substituting 
this in Eq. (Q one finds that different nonzero Fourier modes decouple from each other. This enables an exact 
calculation of any two-point correlation function. For example, as shown in the Appendix, one finds for any k ^ 0, 
(h(k, t)h(k', t)) = [D' /TLk 2 }Sk+k',o in the stationary limit t — > oo. Note that the k = mode is identically zero at 

all times, h(0,t) = 0, which follows from the sum rule jj^ h(x, t)dx = as the height h(x, t), by definition, is always 
measured with respect to its spatial average. Since Eq. yj is linear, the height field h(x, t) is Gaussian for all x and all 
t. Using the result for the two-point correlator mentioned above, one can then easily write down the joint probability 
distribution of the Fourier modes in the stationary state, 



P 



{h(k)} 



oc exp 



T 1 T 

-—*£k 2 h(k)h(-k) 



S[h(0) 



(8) 



where the delta function on the right-hand side of Eq. (JSJ takes care of the "zero-area" constraint. In terms of the 
actual height field h(x,t), the stationary joint distribution becomes [TsL ITflj 



P a t [{h(x)}] = 2V^DL 3/2 exp 



■— /V- 

4DJ \dx 



S [h(0) - h{L)] 5 / h{x)dx 



(9) 



where D = D' /2T and the normalization constant 2\J nDL 3 ^ 2 , ensuring that the joint distribution is normalized, can 
be calculated explicitly |19j. Two delta functions in Eq. © take care respectively of the periodic boundary condition 
h(0) = h(L) and the zero-area constraint. The stationary height distribution at any fixed point x in space is, by 
translational invariance, independent of x and is also a Gaussian 



Pst{h) 



,-h 2 /2(h 2 ) 



(10) 



7 



where the variance (h 2 ) = D'L/12T = DL/6 can be easily computed (see the Appendix). 

Note that in the standard literature on interfaces, one often ignores the "zero area" constraint in the stationary 
measure |2£j- This is justified if one is interested in calculating physical properties in an infinite (L — > oo) system 
where the zero mode h(k = 0) does not play any important role. Besides, in the calculation of certain observables 
even in a finite system, such as the average width in the stationary state or the distribution of the square of the 
spatially averaged width [2lj. the k = mode drops out of the calculation. However, as pointed out in Refs. [l8LIT^ |. 
the "zero area" constraint certainly plays a very crucial role in the calculation of, for example, the distribution of the 
maximal height of the interfaces in a finite system. We will see below that the "zero area" constraint does indeed 
play an important role also in the calculation of spatial survival probabilities in a finite system. 

From the expression of the stationary measure in Eq. © it is evident that stationary paths are locally Brownian, 
i.e., evolve in space as, dh{x)/dx = £(x), where £(x) is a Gaussian white noise with zero mean and a correlator, 
(£(x)£(x')) = 2DS(x — x'). For the periodic boundary condition, the stationary path in space is, in fact, a Brownian 
bridge over x 6 [0, L] that starts at ho = h(0) at x — and ends up at the same point h{L) — ho at x = L. In addition, 
this Brownian bridge has a total 'zero area' under it. It turns out to be convenient to perform the calculations using 
the standard notations of a 'temporal' Brownian motion, dx/dt = £(t) with (£(t)) = and = 2DS(t — t'). 

At the end of the calculations, one can translate back the results to the interface problem upon identifying (i) the 
height of the interface with the position of the 'temporal' Brownian motion, i.e., h = x (ii) the space in the interface 
problem with the time in the 'temporal' Brownian motion, i.e., x = t. Thus, in this notation, the 'temporal' Brownian 
bridge starts at the initial position xo (= ho) at t = and ends at the same position x after a time interval t = T 
(T = L), and enclosing under it a total 'zero area', i.e., J Q T x(t)dt — 0. 

With these notations set up, we first describe a calculation of the spatial survival probability in which the "zero- 
area" constraint is not taken into account. Although the survival probability can be calculated exactly in this case, 
the resulting analytic expression does not show good agreement with numerical results, implying that the "zero-area" 
constraint is crucial for a correct description of the first-passage properties. We then show that the "zero-area" 
constraint can be taken into account in a formally exact path integral treatment. This treatment, however, does not 
lead to a simple closed- form expression for the survival probability that can be compared with numerical results. We 
then use a "deterministic" approximation to obtain closed-form expressions for the survival probabilities and show 
that the analytic results obtained this way provide a reasonably good account of the numerical results. Note that 
since the stationary measure of the Id KPZ interface for the periodic boundary condition is the same as that in the 
EW interface]^, all our steady state results will be valid for the Id KPZ interface as well. 



A. Survival Probability without the Zero-area Constraint 



Let us first recall some basic results for the ordinary free 'temporal' Brownian motion. Consider a Brownian motion, 



§ = «<>• <"> 

where is a zero mean Gaussian white noise with correlator (£(t)£(i')) — 2DS(t — t'). The free propagator of the 
Brownian motion G(x, t\xo, to) defined as the probability that the particle will reach x at time t starting from Xo at 
to can be easily obtained by solving the Fokker-Planck equation, 

d t G(x, t\x ,t ) = Dd 2 x G{x, t|x ,t ), (12) 

subject to the initial condition, G(x, to\xo, to) = S(x — Xo) and the boundary conditions that G — > as x — > ±oo. The 
well known solution is given by, 

G(x,t|x ,io) = , ) e -(-*°) a /*P(«-fa), (13) 

y/4TTD(t - t ) 

valid for all t and to and x and xq- We now ask: what is the probability that the particle reaches x at time t, starting 
at xo at to, but without having crossed the zero in between? This probability P(x,t\xo,to) can be easily calculated 
by solving the same Fokker-Planck equation, but now adding an absorbing boundary condition at x = 0, i.e., insisting 
that P(0, t\xo,to) = for all i|22(. The solution can be easily obtained by the image method and is given by, 

P(X, t\x ,t ) = 1 L-(*-*o) 2 /4C(t-*o) _ e -( X+X0 ) 2 /4D(t-t )] (w) 

y/4ivD(t - t ) L J 



8 



Evidently, this solution satisfies the absorbing boundary condition at x = 0. 

Now, let us consider a Brownian bridge over the interval [0, T]. This means a Brownian motion that starts at xq 
at t = and ends up at the same point xq at time T. Let us ask: what is the probability Q(t, T|a;o) that this process 
(conditioned to be at xo at the two endpoints) does not cross zero in the interval [0, t] where < t < T? To calculate 
this probability, let us divide a typical path of the process over two intervals: [0, t] and [t, T]. Over the first interval 
[0,i], a typical path starts at xo at the left end of the interval and lands up, say at x (where x is a variable) at i, 
without having crossed the zero over the interval [0, t]. This probability for the left interval is (using Eq. 



P L (x,t\x o ,0) 



1 



-{x-xof/iDt _ ~(x+x ) 2 /4Dt 



(15) 



Over the second interval [t,T], the path starting at x at t (left end of the interval [t, T]) reaches at xq at T, but 
there is no restriction over this second interval (the path is allowed to cross zero in this second interval). Thus, this 
probability for the right interval is obtained from the free propagator in Eq. I|13l) . 



P R {x ,T\x,t) 



y/4nD(T-t) 



~{x -x) 2 /4D{T-t) 



(16) 



Due to the Markovian property of the walk, the left and the right intervals are independent. Hence the net probability 
is just the product of the two probabilities, integrated over the position x at the intermediate point t over x € [0, oo]. 
This will be the total probability that a path starting at xq at t — will end up at xq at time T, without having 
crossed the zero in the interval [0, t]. To get the conditional probability Q(t, T\xq) (conditioned that the two ends are 
already given to be at Xo), we need to divide this probability by the factor 1/ 'v 'AttDT which is just the probability of 
a free path landing up at xq at T, starting at xq at time 0. Thus, finally, we get 



Q(t,T\x ) 



T 



dx 



ATtDt{T - t) J 
This integral can be done in closed form and one gets, 



e -(x-x y/4Dt _ e -(x+x y /ADt 



Q(t,T\x ) = i 



1 + erf xq 




T 



ADt{T - t) 



e-*S/ OT erfc 



xo 



-(x -xf/4,D(T-t) 



(T - 2t) 



y/ADtT(T - t) 



(17) 



(18) 



where erf (a;) = -4= e u du is the error function and erfc(x) = 1 — erf(x). 

So, now we can interpret these results in terms of the stationary state of the EW interfaces. Identifying xq = ho 
and T = L and t = x, Q(x, L\ho) in Eq. I|18[l is just the probability that the stationary interface, given its height 
at the two ends of the sample, does not cross zero in the spatial interval [0, x] and is given by, for ho > 0, 



Q{x,L\h ) 



1 



1 + erf h, 




e -ht/DL elic ho _ 



(L — 2 x) 



ADx(L-x)J ~ y/4DxL(L-x) 

Now, we need to average Q{x, L\ho) over the stationary distributions of ho given in Eq. (|10|) . 

Q(x, L) = 2 dho Q{x, L\ho) P st (h ). 
Jo 



(19) 



(20) 



The factor 2 comes from the fact that the stationary distribution becomes twice its value if restricted over only the 
positive half-space ho € [0, oo]. The integral in Eq. (|20|l can be done in closed form. We need to use the identity, 



crf(zx) = —tan 1 (z), 



(21) 



which can be easily proved by differentiating both sides with respect to z, performing the resulting integral and then 
integrating back with respect to z. Our final result is: Q(x,L) = f(x/L) for all x and L, where the scaling function 
f(u) is given exactly by, 
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FIG. 5: (Color online) Comparison of the analytic result for the spatial survival probability (Eq. 1221 with u = x/L, Q(x,L) = 
f(u)), obtained without enforcing the zero-area constraint ((red) dashed line) with the numerical result for L = 10 4 ((black) 
solid line). 



One can easily check that /(0) = 1 as it should be. Also, note that /(l) = 1 — V~3/2 is nonzero. Interestingly, this 
scaling function f(u) does not depend on the system parameters D' or T. 

The function f(u) in Eq. 121'P is plotted vs. u in FigEJ and compared with the numerical result obtained for 
L = 10 4 . The agreement between the analytic and numerical results is not satisfactory. In particular, the numerical 
result for /(it) goes to very small values as u is increased above 0.5, while the analytic curve shows a finite value 
f(l) = 1 — v3/2 even at u = 1. It is, therefore, clear that the "zero-area" constraint has to be included in the 
calculation for a correct description of the survival probability. 



B. Survival Probability with the "zero-area" Constraint 



In the previous section, we did not take into account the constraint that the total area under the Brownian bridge 
x(t) going from xq a t time r = to Xo at time r = T is actually zero. In this subsection, we perform the calculation 
taking into account this "zero-area" constraint. 

We first define Qo(t,T\xo) to be the probability that the process x{t) starting at xq at r = does not cross zero 
up to time t (see Fig. EJ) where < t < T, given that the process ends up at xo at t — T and that the total area under 
the process from to T is 0. The subscript in Qo indicates that the fact that we are considering only those paths 
(out of all possible paths starting at xq) whose area is 0. Note that the analogous quantity Q(t, T\xq) without the 
area constraint was computed in Eq. I|18|) . Formally, one can express this probability Q (t,T\xo) as follows, 

£*H) (23) 

where the angular brackets (• ■ • ) denote an average over all possible Brownian paths that start at xq at time 0. The 
numerator Af in Eq. I|23|l represents the joint probability of 3 events: (i) the probability that the path does not cross 
zero upto t (the first factor) (ii) the probability that the path finally ends up at xq at time r = T (the second factor) 
and (iii) the probability that the area under the process upto time T is zero (the third factor). The denominator 
V in Eq. i|23*|) represents the joint probability of events (ii) and (iii). So, the ratio J\f/T> represents the conditional 
probability Qo(t, T\xq), i.e., the fraction of paths that start at xq and satisfy (i), (ii) and (iii), out of all paths that 
start at Xg and satisfy (ii) and (iii). 

This numerator can further be split into two parts: (a) the probability that a path starting at Xo reaches x at the 
intermediate time t without crossing zero and enclosing an area, say A > (the area is positive since the path does 



'(s(t)) 6[x(T)-x }6 



(s \x(T) — xq] S J Q T x(t) dr 
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FIG. 6: A typical Brownian path starting at xo and reaching x at time t without crossing the origin and then freely propagating 
from x at t to xq at time T. 



not cross zero in [0,t]) and (b) the subsequent probability that the path starting at x at time t reaches xq at time T 
and enclosing an area —A. This ensures that the total area upto T is zero. However, we then have to integrate over 
all possible values of A > and x > 0. Thus, one can rewrite Eq. 1|23[) as, 



] (t,T\x ) 



J °° dA J °° dx G+(x, A, t\x ,0, 0)G(x , -A, T - t\x, 0, 0) 
G(x o ,0,T|a:o,0,0) 



(24) 



where we define G(x, A, t\xo, Aq, 0) to be the probability that the joint i two-variable Gaussian process [x(t),A(t) = 
Aq + J x(r)dr)] (i.e., the Brownian curve itself and the integral under the Brownian curve) reaches [x,A] at time t, 
starting from the initial value [xq, Aq] at time 0. So, G(x, A, t\xo, Aq, 0) is just the propagator for this joint Gaussian 
process. Note that this process is free in the sense that it has no restriction for x(t) to be only positive. Clearly, the 
denominator V in Eq. I)23|l is just G{xq, 0, T\xq, 0, 0) since, by definition, this quantity represents the probability the 
process [x(t), A(t)] will reach its final value [xq, 0] (note that the final area is zero) at time T, starting from its initial 
value [ieoiO] (since the initial area is zero). Similarly, the part (b) in the numerator is just G(xq, — A, T — i|x,0, 0) 
since that represents the joint probability that a path starting at x at time t and initial area will end at xq at time T 
with area — A (note that we have made a shift [t, T] — > [0, T — t] which is allowed due to time translational invariance 
of Brownian motion). Finally we define G + (x, A, t\xo, 0, 0) to be the joint probability that the process [x(t), A(t)] will 
reach [x, A] at time t starting from [acojO] and without x crossing the origin in the interval [0, t]. 

Now, the free propagator G(x, A, t\xo, Aq, 0) is easy to compute. This is just the joint bivariate Gaussian distribution 
of the random variable x(t) and A(t) = Aq + J Q x(t')dt' . Their correlation matrix can be easily computed and the 
propagator is just the exponential of the inverse correlator. Indeed, this result is already well known, since this is 
just a random acceleration problem: dA/dt = xU) and dx/dt = rj(t) which implies d 2 A/dt 2 — rj(t). The result for the 
propagator can be found, for example, in Ref . fl3. l23| . We just quote this result here, 



y/3 

G(x, A, t\x ,A Q , 0) = 27rDi _ 2 exp 



1 f 3 1 

— i^(A-A - xt){A -Aq- x t) + -{x - Xq) 2 



(25) 



Substituting this exact propagator in Eq. I|24|l . we get 

2 



h{t,T\x ) = 



T 



T-t 
x exp 



6 3x /dt dA dxG + (x,A,t\x ,0,0)x 
Jo Jo 

(A + x(T - *)) (A + x (T - t)) - * (a; - x f 



D{T-tf 



D(T - t) 



(26) 
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Note that we still have to determine the restricted propagator G + (x, A, t\xo, 0, 0) and then we have to average 
Q(t, T\xq) over the stationary distribution of Xq to calculate Qo{t, T) = 2 J Q dxoQo(t, T, xo)P st {xo), where the distri- 
bution Pst(xo) is given in Eq. ijlH)!. With the identification ho = x Q and L = T, one gets 



PsM = \l^e-^/ DT , (27) 



which gives 



Q (t,T) = 2\l^ I dx e- 3x o/ DT Q Q (t,T,x ). (28) 



Substituting the expression of Qo(t, T\xq) from Eq. (|2lj|) in Eq. we find that the factor e 3x l/ DT cancels out and 
we get a formally exact result, 



x exp 



(A + x{T - <)) (A + x (T - t)) - - x f 



D(T-tf y v //v v " D(T — t) 



(29) 



So the remaining task is to evaluate the restricted propagator G + (x, A, t\xo, 0, 0) which we do in the next two sub- 
sections. 



Exact Calculation of the Restricted Propagator 



We note that the restricted propagator G + (x, A, t\xo, 0, 0) is just the joint probability that the process x(t) starting 
at xq at time reaches x at time t without crossing the origin and that the area under the curve is A. Thus it is 
given by, 



G + (x,A,t\x , 0,0) 



II<?(x(r)) 



T=0 



5 [x{t) 



x(t) dr — A 



(30) 



where the angular brackets (• • • ) denote an average over all possible paths starting at xo at time 0. Note that if 
x = x = 0, this restricted Brownian process is just a Brownian excursion over the interval [0, t] and the propagator 
G + (0, A, t\0, 0, 0) is just the (unnormalizcd) probability density of the area under a Brownian excursion. A Brownian 
excursion is simply a Brownian path that propagates from x(0) — to x(t) = over [0,t] but is conditioned to stay 
positive in between. The probability distribution of the area under a Brownian excursion was calculated exactly and 
is known as the Airy distribution function, a complicated function that involves the zeros of the Airy function but is 
not the Airy function itself |24l 123 . l2o. f27L l28j . This Airy distribution function has been a subject of intense study for 
the past several years as it has resurfaced in many problems in computer science p^. l27| . graph theory 29], and two 
dimensional polygon problems |30j|. Recently it was shown that the same Airy distribution function also describes the 
distribution of maximal height in the stationary state of fluctuating interfaces [l8l ITsj. In Ref. 0]> a path integral 
derivation of the area distribution under a Brownian excursion was provided. In our present problem, we have a 
generalization of this problem where a Brownian path propagates from xq at r — to x at r = t, staying positive 
in between. The nonzero values of the initial and the final positions, x$ and x, make explicit calculation of the 
distribution of the area under such a path difficult, as demonstrated below. Here we follow the path integral method 
used in Ref. generalized to nonzero xq and x. 

Let us define the Laplace transform of G + with respect to the area A, i.e. 



G(x,x ,t,\) 



G + (x, A,t\x , 0,0)e 



— XA 



dA. 



(31) 



Note that since the path is restricted to be on the positive side, the area is always positive and hence a Laplace 
transform (rather than a Fourier transform) is more suitable. Taking Laplace transform of Eq. (|30|l we get 



G(x,x ,t, A) 



i[0(x(r)) 



5 [x(t) 



x e 



" A Jo' x i T ) dT 



(32) 



12 



Using the Brownian measure of the paths in Eq. G(x, xq, t, A) can be expressed as a path integral, 



G{x, Xq, t, A) oc 



x{t)—x 



T>x(t) 9[x(t)] exp 



dr 



1 (dx{r) 



AD \ dr 



Xx(t) 



(33) 



' x(0)=x o 

Using the bra-ket notation, one can re-express the right hand side of Eq. (|33|l as a quantum mechanical propagator, 

G(x,Xo,t,\) = (x Q \e~'^ t \x}, (34) 
where the Hamiltonian H can be written (in the position basis), 

d 2 



H = -D^ + V(x), 



(35) 



where the quantum potential V(x) — \x for x > and in addition V(0) = oo (i.e., there is a hard wall at the origin) 
which takes into account the fact that all paths in Eq. Q33[) are restricted to be on the positive side and can not enter 
the region x < 0. Expanding the right hand side of Eq. (|34|) in the eigenbasis of the Hamiltonian H we get, 



x e 



-El 



G(x,x ,t,X) = Y,{x \E)(E\s 

E 

= ^2i> E (x)r E ^o)e- Et , 



(36) 



where E denotes the eigenvalues of H and the eigenfunction ipE(x) satisfies the Schrodinger equation in the region 

x € [0, oo], 



d 2 

-D—t/j(x) + Xxip(x) = Eip(x), 
ax z 



(37) 



with the boundary conditions, ^>(oo) = and i{}(0) = 0, the latter reflecting the hard wall at the origin. Making the 
change of variable, z = {\/ D) 1 / 3 {x — E/X), one can recast the Schrodinger equation as the following Airy differential 
equation, 



d 2 ip 
Ih 2 



zi) = o, 



whose general solution is given in terms of two Airy functions, 

ip(z) =£Ai(z) + CBi(z), 



(38) 



(39) 



where B and C are arbitrary constants. The function Bi(^) diverges as z — > oo (see Ref . |31| ) indicating C = 0. Going 
back to the original x variable, we have 



ip(x) = B Ai 



1/3 



The other boundary condition ip(x — 0) = determines the eigenvalues, 

E 



Ai 



(DX 2 ) 1 / 3 



0. 



(40) 



(41) 



It is known (see Ref. [3ll] . page 446) that the Airy function Ai(x) has zeroes on the negative x axis at x = —ai. For 
example, «i = 2.33, a 2 — 4.08, a% — 5.52, — 6.78 etc. Thus we get the exact eigenvalues 



E^a^DX 2 ) 1 ^. 

The amplitude B in Eq. I|40|l is determined from the normalization, J °° \ipE(x)\ 2 dx = 1 and we get 

1 



B,V = 



j^Ai 2 [(X/D^y-at] dy' 



(42) 



(43) 
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Finally, putting everything back in Eq. (|36|) we get the exact Laplace transform of the restricted propagator, 

Aip/D) 1 / 3 ^ - a,} Ai[(A/Z?) 1 /3 x _ at 



G(x, x ,t, A) = ^2 ' 



-aiD L 



/ °°Ai 2 [(\/Dy/3y-a t ]dy 



(44) 



where — o^'s are the zeros of the Airy function Ai(x). Formally inverting the Laplace transform using the Bromwitch 
formula, we get the exact expression for the restricted propagator, 



G + (x,A,t\x o ,0,0) 



Ao+ioo ^\ 

— e XA G(x,x ,t,X), 

Ao-ioo /7ri 



(45) 



where G{x, xq, t, A) is given by Eq. Q44JI and the integration is along any imaginary axis whose real part Ao must be 
to the right of all singularities of the integrand. 

Substituting this restricted propagator from Eqs. (|45|l and (|44|l in Eq. I)29|) gives us the formal exact answer of the 
no zero crossing probability Qo(t, T) which, one can easily check, has the scaling form Qo(t, T) = f(t/T). However this 
formal expression of the scaling function f(u), though exact, is hardly useful to make comparison with the numerical 
data. This is because, we can not invert the Laplace transform explicitly in Eq. 14511 (even after that one needs to do 
the 3 integrals over xo, x and A in Eq. i|29|) which looks hopeless!). 

So, the question is: Can one find an approximate way to estimate the restricted propagator G + (x, A, t\xo, 0, 0) 
and then use this approximate result in Eq. (j29(l . do the 3 integrals and derive an explicit expression for the scaling 
function f(u)7 In the next subsection, we indeed perform these steps and derive an explicit expression of f(u) which, 
though not exact, is expected to be quite good. 



2. A Deterministic Approximation for the Restricted Propagator G 
As a first step, let us rewrite the restricted propagator G + (x, A, t\xo, 0, 0) in Eq. as follows, 



G + (x,A,t\x , 0,0) = 



< 




8 [x(t) -x]8 


Jo x(t) dr- A 


) 


< 


l£=o" (*(*-))_ 


S[x(t) -x}\ 



X\9(x{t)) 



T=0 



5 [x(t) — x] 



= W[A,x,x o ,t}P(x,t\x ,0), (46) 

where we just divided and multiplied the right hand side of Eq. I|30|) by the same factor P(x, t\xo, 0) — 
^ Ilt=o ^ ( X ( T )) $ l x (t) ~~ x ]\ This factor is simply the probability that the path starting at xq at time reaches x 
at time t without crossing the origin in the interval [0, t] and has already been calculated by the image method in 
Eq. lO. Thus, we have, 



P(a;,t|a:o,0) 




5 [x(t) - x] 



) 2 /iDt _ e -(x+x„) 2 /4Dt 



(47) 



The quantity W(A,x,x ,t) represents the expression inside the first parenthesis on the right hand side of Eq. l|4^|) 
which we can write as 



W(A,x,x ,t) 





Jo x( T ) dr- A 






5[x(t)-x]\ 


( 




(x(r))_ 


5 [x(t) ~ 


x)) 



(48) 



which is simply the probability distribution of the the area under the process in [0,t], given that the process reaches 
from xq to x in time t without crossing the origin in between. 

So, if we can estimate this conditional area distribution W(A,x,xo,t), then knowing the exact P(x,t\xo,0) from 
Eq. l(4T|) . we will be able to estimate the restricted propagator G + (x, A, t\xo, 0, 0) from Eq. (J3BJ. 

Note that so far we have not made any approximation. To estimate the conditional area distribution PF(j4, x, Xq, t) 
defined in Eq. Q4H|). we now make a "deterministic" approximation as follows. Note that W(A,x,Xo,t) is just the 
fraction of paths that have an area A, amongst all possible paths that go from xq to x without crossing the origin in 
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between. Now, the paths that go from xq at r = to x at r = t without crossing the origin in between are likely to 
cluster around an optimal path, i.e. most of these paths lie "close" to an optimal path. This optimal path is the one 
which has the highest probability amongst all possible paths going from xq at r = to x at r = t without crossing 
the origin in between. Assuming the existence of such an optimal path x op t(r), the "deterministic" approximation 
consists in writing 

W(A, x , X) t)Ks(A-J x opt (r)dA . (49) 

Thus, within this aproximation we ignore the fluctuations that arise from non-optimal paths. 

The next step is to actually find the optimal path Xopt(r), i.e. the path with the highest probability, amongst all 
possible paths that satisfy the following constraints: (i) x(0) = xq; (ii) x(t) — x and (iii) x(t) > for all < r < t, i.e. 
the path stays positive in the interval r £ [0, i]. One knows from the principle of least actions that the optimal path 
is the so called "classical" path that satisfies Newton's equation of motion. In our problem, the action in the path 

integral, S = jjj 0^§p^ j dr, corresponds to that of a free particle. So, the optimal path satisfies the Newton's 

law for a free partricle, namely d 2 x/dT 2 = 0, starting from x(0) = xq and ending at x(t) = x for all < r < t. The 
solution is trivially, 

£opt(V) == x + (x - Xq)~- (50) 

Note that this solution automatically satisfies the condition (iii) mentioned above, i.e. it stays positive in the interval 
r G [0,£]. Now, the area under the optimal path is simply 

f* 1 

^opt = J X opt (r)dT = - [Xq + x\t, (51) 

Before proceeding to the calculation of the survival probability with this optimal choice, it is instructive to ask 
how the calculated survival probability depends on the choice of the deterministic path. In other words, within the 
"deterministic" approximation in Eq. I|49l) , how does the final result vary if instead of the optimal path in Eq. I|5U[) , we 
use some other paths? To test this, we actually consider a one parameter family of paths that satisfy the constraints 
(i), (ii) and (iii) above. This family of paths is characterized by a single parameter (3 > 0, 

X/3 (t) = X + [X - Xq) (-J . (52) 

Clearly, as shown above, (3 = 1 corresponds to the optimal path. In the following, we will however calculate the 
survival probability for all j3 > to test how sensitive the final answer is to the optiaml choice (3—1. We will see, 
somewhat surprisingly, that the final scaling function /(it) characterising the survival probability depends only very 
weakly on [3. 

The area under the deterministic path in Eq. I|52|) is simply 

[ l (3 

A((3, Xq, x,t) — / xp{r)dT — [axQ + (1 — a) x]t ; where a= -. (53) 

Jo 1 + (3 

The optimal path corresponds to the choice (3—1, i.e. a — 1/2. Our approximation in Eq. I|49() then reads 

W(A,x,x ,t) w 5 [A - (ax + (1 - a) x)t] , (54) 

the optimal choice being a = 1/2. 

Within this deterministic approximation, we then have an explicit expression for the restricted propagator [on 
substituting the results in Eqs. l|T7|) and (|51Jl in Eq. I0BJ|], 

G + (x,A,t\xQ,0,0) « S lA-(axo + (l-a)x)t] l_ (x _ Xof/ , m _ e - {x+Xo?/iDt ^ (K) 



3. Explicit Expression for the Scaling Function f(u) using the Deterministic Approximation 



On substituting the expression of G + (x, A, t\xo, 0, 0) from Eq. I)55|l into Eq. (|29|) . we can do the integral over A 
trivially, since it involves a delta function. Inside the exponential on the right hand side of Eq. I|29|l , we have to replace 
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A by (axo + (1 — a)x)t. This gives, 
Qo(t,T) = W3a/1 



(£0 


2 


f°° dx 


f°° dx f 




y/AnDt J 





e -(x-x y/±Dt _ e -(x+x y /ADt 



x exp 



D(T-t) 1 



(a(xQ — x) t + xT) ((1 — a) (a: — icq) t + T) 



D(T - i) 



(x - x Q y 



(56) 



We next define the scaling variables: y = xo /V ADt, z — x/yADt and u = t/T. Then Qo(t,T) in Eq. (|56|) becomes 
only a function of it = t/T, i.e., Qo(t,T) = f(t/T) where the scaling function /(it) is given from Eq. ifSBl . 



/(«) = 



4V3u 



OO /*OG 



tt(I-w) 2 J 
x exp 



dy / dz 



e -(v-*Y - e -(y+ z r 



Au 



(i-uy 



{ (7U 2 + (3a - 2)u + 1) y 2 + (ju 2 + (1 - 3a)u + l) z 2 + (l + u - 2 7 u 2 ) yz} 



where 7 = 1 — 3a + 3a 2 . The right hand side of Eq. (|57|l can be reorganized as, 



where 



h{u) 
h(u) 



/(«) 



OO /*OG 



4V3u 
tt(1 — u) : 



[h(u)-I 2 (u)}. 



(57) 



(58) 



Jo 

OO /"OO 



dydz exp [— r(u)y 2 — s(u)z 2 +p(u)yz\ , 
dydz exp [— r{u)y 2 — s(u)z 2 + q(u)yz] , 



(59) 



'0 Jo 

where r(u) = (3(l-4a + 4a 2 )u 3 + (12a-5)w 2 + u+l)/(l - u) 3 , s(u) = (3(l-4a + 4a 2 )u 3 + (7-12a)u 2 + u + l)/(l - u) 3 , 
p(u) = 2(3(1 - 4a + 4a 2 )u 3 + u 2 - 5it+ 1)/(1 - u) 3 and q(u) = -2((-5 + 12a - 12a 2 )w 3 + 5w 2 - u + 1)/(1 - u) 3 . To 
do the double integrals in Eq. I|59(l . it is convenient to scale y = Y/yfr{u) and z = Z/y/s(u). This gives, 

(1-n) 3 ^ ^ 



7 x (u) = 
M = 



\/B(u) Jo Jo 

a\3 poo poo 
— U)° ' ' 



dYdZexp [-(Y 2 + Z 2 + 2A 1 {u)YZ] 

poo poo 

/ / dYdZ exp [- (Y 2 + Z 2 + 2A 2 [u)YZ] 
Jo Jo 



(60) 



where 



B(u) = (3(1 - 4a + 4a 2 )u 3 + (12a - 5)w 2 + u + 1)(3(1 - 4a + 4a 2 )u 3 + (7 - 12a)u 2 +u+l) 
3(1 - 4a + 4a 2 )u 3 + u 2 - 5u + 1 



Ax(u) 
A 2 (u) 



(-5 + 12a - 12a 2 )u 3 + 5u 2 - u + 1 



We next use the identity, 



/ / dYdZcxp[-{Y 2 + Z 2 + 2AYZ)] = - - 
Jo Jo 4 7o 1 



dB 



+ Asin(0) ' 



(61) 



(62) 



which can be easily established by going to the polar coordinates Y — Rcos(8) and Z = i?sin(6*) and by performing 
the radial integration. The integral J(A) = f£ 1+ Asi n (e) can ^ e done in closed form and one gets, 



J(A) = 



1 



VT^A 2 
1 



[7r-2sin~ 1 (yl)], for |A| < 1, 
log " ' y - ^= = for|A|>l. 



A - VA 2 - 1 



(63) 
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FIG. 7: (Color online) Comparison of the analytic result for the spatial survival probability (Eq. I)(i4^ with u = x/L, Q(x,L) = 
f(u)), obtained with the zero-area constraint and using the deterministic approximation with a — 1/2 ((red) dashed line) with 
the numerical result for L — 10 4 ((black) solid line). The analytic result for f(u) obtained using a — 0.1 is also shown ((blue) 
dotted line) for comparison. The inset shows the same three plots for larger values of u, using a logarithmic y-scale to illustrate 
the dependence of f(u) on the parameter a for relatively large values of u. 



Putting all these results back in Eq. (JSSJl, we get an explicit result for the scaling function f(u), 

= V3 ^L-u) [J{Mu)) - J{ A 2 (u))} , (64) 

where the functions J, A\, A2 and B are given respectively in Eqs. Ij63(l and H61(l . This is our main result, obtained 
using the deterministic approximation where we keep only the contribution from the optimal path but ignore the 
fluctuations around the optimal path. 

The function in Eq. (|64|l can be easily evaluated for all < a < 1. Note that, by definition, the a dependence of 
f(u) is symmetric about a = 1/2, so it suffices to use < a < 1/2 with a = 1/2 being the optimal choice. In Fig. [7| 
we compare the analytically obtained f(u) corresponding to the optimal choice a = 1/2 with the actual f(u) obtained 
via numerical simulation on a latice of L = 10 4 sites. The analytical scaling function f(u) seems to compare well with 
the numerical one, given especially the fact that there was no fitting parameter involved. 

The function f(u) turns out to be rather insensitive to the value of the parameter a. In Fig. \7\ we have also plotted 
f(u) for the choice a = 0.1. The f(u) obtained for this non-optimal choice of a is virtually indistinguishable from 
that obtained for the optimal choice, a = 1/2. This can be understood from the asymptotic behavior of /(it) near 
u = 0. As u — * 0, one can show that, 

f{u) = ! _ ±/3« + 8 m 3 /2 + 473 u5/2 _ 32V3ag~a) u7/2 + 9/2 

7T 7TV3 7T 7T 

Note that the a-dependence appears only in the 5th term in the small-u expansion, showing that the function f(u) 
is highly insensitive to a for small u. Since f(u) decreases rapidly with increasing u, the dependence of f(u) on a for 
relatively large values of u is not visible in the linear plot in Fig. This dependence is evident from the plots in the 
inset of Fig. [7| where the results for f(u) for two different values of a (the optimal value, a = 1/2 and a = 0.1) arc 
shown in a logarithmic scale and compared with the numerical data for L — 10 4 . 

The dependence of f(u) on the parameter a for relatively large values of u is clearly seen in the asymptotic behavior 
of f(u) as u — > 1. One finds that in powers of e = 1 — u where e — > 0, 

f(u = 1 - e) = c 4 (a)e 4 + c 5 (a)e 5 + c 6 (a)e 6 + 0(e 7 ), (66) 
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where the coefficients 04(a), 05(a), c§{a) etc. are complicated functions of the parameter a. In particular, for the 
optimal case a = 1/2, we get 

/(«=!-£) = — Ue 4 + ^e 5 + —^e 6 + 0(e 7 ). (67) 
9\/37r 3V3tt 270V3tt 



C. Survival Probability for Finite Initial Conditions 

Analytic results for the FIC spatial survival probability Q pic discussed in section [I] may be obtained from the 
calculations described above in the following way. We consider the probability that, given the constraint that the 
height ho at x = lies between — w and +w with w -C W(L), the steady-state width of the interface of size L, the 
interface height does not cross zero within distance x. This probability is normalized by the probability of finding the 
initial height between —w and +w [i.e. by 2 f Q P st (h)dh where P s t(h) is the Gaussian probability distribution for the 
height in the steady state (see Eq. J3J]. In the first calculation without the zero-area constraint, we use Eq. (|19fl and 
integrate the quantity Q(x, L\ho)P a t{ho) over h Q between and w (with a factor 2 to take care of negative values), 
assuming that w/W(L) <C 1, and u — x/L is of order unity. Keeping terms of the lowest (linear) order in y = w/W(L) 
in the expansions of the exponentials and error functions (the latter is justified as long as u is not very close to 
or 1), we get the following result for the normalized probability that the height does not cross zero within distance 
x = uL, when the initial height lies between — w and +w with w = yW(L): 



Qfic{x, L, w) = — u)/u. (68) 

v 247r 

This is clearly consistent with the scaling form of Eq. (0 - the sampling interval 5x is zero in the present continuum 
description. For small it, the scaling function shows a power-law decay with exponent 1/2, fFic{u,0,y) ~ A(y)/y/u, 
with A(y) = y/\f2An. 

In the calculation with the zero-area constraint, we use the results obtained above using the deterministic approxi- 
mation. Specifically, we consider Eq. I|56|) (with t replaced by x and T replaced by L, so that u = t/T = x/L), and 
do the Xq integration between and w instead of between and 00, and then divide the result by [2 J Q P st (h)dh] for 
normalization. Again, keeping terms to the lowest order in y = w/W(L), we obtain the result 

Q FIC ( X ,L, W ) = f FIC (u,0,y) = VTTW^ 1 "^^^"^^, (69) 

where the functions A±(u), A2(u) and B(u) are defined in Eq. (|61|) and r(u) is defined in the line after Eq. I|59|l . We 
use the optimal value, a = 1/2, in evaluating these functions. In the small-u limit, this expression reduces to the 
same form, A(jj)/^/u, as that found in the calculation without the zero-area constraint. This is expected, since the 
zero-area constraint becomes important for values of x comparable to L. 

As shown in Fig.0wb.ere we have plotted the scaling functions according to Eq. (|68[) and Eq. 169fl for w/W(L) = 0.02 
versus u — x/L on a log- log scale, the results obtained with and without the zero-area constraint agree with each 
other for small u, but show differences as it increases above about 0.01. In the same Figure, we have also shown 
the numerical results obtained for the same value of y and L = 10 4 . The numerical data show the expected A/^fu 
behavior for small u, but the coefficient A obtained from a fit to the numerical data turns out to be substantially 
larger than the value 0.02/\/247r ~ 0.0023 predicted by the analytic calculation - the value of A obtained from the fit 
is close to 0.0073. In FiglHI we have also shown the result of Eq. (|69l) multiplied by 3.167, to take into account this 
difference between the two values of A. It is clear from the plots that the analytic result multiplied by this empirical 
factor provides a good description of the numerical data - the level of agreement is roughly similar to that found in 
FigEI for the steady-state survival probability. 

The reason for the necessity of multiplying the analytic result by a numerical factor in order to obtain approximate 
agreement with the numerical result lies in the use of discrete sampling in the numerical calculation of the FIC survival 
probability. In Eqs. I|68() and l|69(l above, the survival probability goes to zero as y — > 0. This reflects the fact that 
in the continuum limit, the probability of not crossing zero over a finite distance x = uL is zero if the initial height 
is zero. This, however, is not true when the sampling of the height is done at discrete points, x n — ndx, where n is 
a positive integer and 8x is the sampling interval, taken to be equal to the spatial discretization scale (5x = 1) for 
most of the numerical results reported here. This is because the probability calculated in the numerical work does 
not take into account the (many) zero crossings that would have taken place between X\ and x-i in the continuum 
limit if h(x\) is very close to 0. Even if the height at the initial point x\ is very close to zero, the probability that 
the height at the next point X2 — x\ + 1 has the same sign as that of the height at the initial point (this is the first 
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FIG. 8: (Color online) Comparison of analytic and numerical results for the FIC spatial survival probability Qfic(x, L, w) with 
w/W(L) — 0.02. The (blue) dashed line and the (red) solid line show, respectively, the analytic results obtained without the 
zero-area constraint (Eq. lEHQ an< ^ with the zero-area constraint (Eq. I|69^ with a = 1/2). The numerical results for L = 10 4 
are shown by (black) circles, and the (magenta) dotted line going through these data points represents the result of Eq. 16911 
multiplied by 3.167. 



value of the measured "no zero crossing probability") is actually close to 0.5 in all the numerical simulations - for 
hi slightly above zero, the probability of h 2 remaining positive is close to 0.5 because, as discussed in section the 
height difference Si = h 2 — hi is a Gaussian random variable with zero mean. 

Our numerical results (such as those shown in Fig^J suggest that for a fixed value of y — w/W(L), the scaling 
functions fFic{u, Sx/L, w) for different values of Sx/L differ from one another mainly by an overall numerical prefactor 
that decreases as the value of Sx/L is reduced. This is why the analytic results for the survival probability show 
reasonable agreement with the numerical data (as shown in Fig|8J when the former are multiplied by a suitable 
numerical factor. An approximate analytic estimate of this numerical factor may be obtained in the following way. 
Using the statistical properties of the height difference variables {sj} mentioned in section [HJ it is easy to show that 
for L ^S> 1, the probability that the height h 2 at lattice site 2 has the same sign as that of hi, the height at site 1, is 
given by 0.5 [1 + erf(/ii/v / 2)]. The "one-step" FIC survival probability for discrete sampling with Sx = 1 is then given 
by 



id(w, L) = Qfic{x = 1, L, Sx — 1, w) = 0.5 



fif dhi exp(-hj/2W)erf(hi/V2) 
J™ dh 1 exp(-h 2 1 /2W) 



(70) 



where W(L) = y/L/12 is the steady-state width of the interface. For w <§; W(L), this becomes a function of w only: 

Q u (w) = 0.5[1 + erf(ti7/V2) - y/2fr(l - e^'^/w]. (71) 

As expected, Qid(w) goes to 0.5 as w — > 0. Our numerical results for Qid{w) are in good agreement with this 
analytic prediction. The "one-step" FIC survival probability in the continuum limit may be obtained by setting 
x = 1 in Eq. I|19l) . integrating Q(x — 1, L\ho)P s t(ho) over ho between and w, and dividing by J Q P s t(ho)dho for 
normalization. For w <C W(L), L 1, this leads to the result 



Qic(w) = Qfic{x = 1, L, Sx = 0, w) = erf(w)/ 



/2/^(l _ e ~ w 2/2 )/w = 2Q ld (w) - 



1. 



(72) 



Thus, Qi c (w) goes to zero as w — * 0, as expected. However, both Qid(w) and Qi c (w) approach unity for L ^> 1, 
w 3> 1, w/W(L) — > 0, indicating that the numerical and analytic results would agree with each other in this limit. 
For w/W(L) = 0.02, L = 10 4 (the values for which numerical results are shown in Fig|S|), the values of the one-step 
survival probabilities are Qid — 0.612 and Qi c = 0.224. The ratio of these two number is 2.73, which is slightly smaller 
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than the empirical multiplicative factor used in Fig|S]to bring the analytic result for the FIC survival probability in 
approximate agreement with the numerical data. This difference reflects the fact that the empirical value used in 
FiglHlwas obtained by considering the analytic and numerical results for a range of values of x, whereas the analytic 
estimate of the multiplicative factor is obtained by considering only one point. 

IV. SUMMARY AND DISCUSSIONS 

In this paper, we have presented analytic and numerical results for the spatial survival probabilities for Id EW 
interfaces in the steady state. In Id the same steady state results also hold for the KPZ interface. We have studied 
both the steady-state and the FIC survival probabilities. The numerical results show that these survival probabilities 
exhibit simple scaling behavior as functions of the system size and the sampling interval used in the measurement. In 
the analytic work, we have used a "deterministic" approximation to obtain closed-form expressions for these scaling 
functions from an exact path integral treatment of a mapping of the problem to Id Brownian motion. The analytic 
results show fairly good agreement with the numerical data without having to use any adjustable parameter. The 
remaining differences between the analytic and numerical results may be attributed to (a) the use of the determin- 
istic approximation in obtaining the analytic results, and (b) the use of a finite sampling interval in the numerical 
calculations. The effect of discrete sampling is small in the case of the steady-state survival probability. For the FIC 
survival probability, the dependence of the numerical results on the value of the sampling interval used in the mea- 
surement is approximately described by an overall multiplicative factor. The value of this multiplicative factor can be 
estimated analytically by considering the one-step survival probability. Further analytic work on the dependence of 
the survival probabilities on the sampling interval would be interesting and useful, especially because measurements 
of these probabilities in simulations and experiments always involve discrete sampling. While some progress in this 
direction has been made ^(| , a complete analysis of the effects of discrete sampling remains a challenging theoretical 
problem. 

On the experimental side, fluctuating steps on a vicinal surface provide a physical realization of Id EW interfaces 
if the kinetics is dominated by attachment /detachmentprocesses 1 1 5| . While experimental studies of temporal persis- 
tence and survival probabilities have been carried out [y, 0, 01 El f° r this system, we are not aware of any experimental 
investigation of spatial first-passage properties of fluctuating steps. Such investigations would be most welcome. In 
other experimental systems such as combustion fronts in paper which also involve Id interfaces described by the EW 
or the KPZ equation, the spatial persistence has been recently investigated ^3- ^ would be interesting to see if our 
theoretical predictions on the spatial survival probability can also be verified experimentally in such systems. 
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APPENDIX A: FINITE ID EW INTERFACES WITH PERIODIC BOUNDARY CONDITION 



Consider the EW equation Q on a finite line of size L with periodic boundary condition, h(x+L, t) = h(x, t)). Since 
h(x,t)) is a periodic function, one can decompose it into a Fourier series, h(x, t) = J2 k H k > t)e lkx where k = 2-KmjL 
with m = 0,±1,±2, .... Thus 



h(x,t)= h(m,t)e 2nl 



c/L 



(Al) 



m— — oo 



where the Fourier coefficients h{m, t) are given by the inversion formula, 



h(m, t) 



L Jo 



h(x,t) e - 27Ilmx/L dx. 



Substituting Eq. (|A1| in the EW equation £Q) one gets, 



- 4n 2 m 2 
Othim, t) = — — T him, t) + Tj(m, i), 



(A2) 



(A3) 



for all m ^ 0, where (rj(m,t)) — and {fi{m,t)fi{m! , i') = =^-S m + m i ,o with S m . n being the Kronecker delta function. 
Note that for m = 0, h(0,t) = at all i, because of the sum rule, J Q L h(x,t)dx — 0. Solving Eq. I|A3(1 with initial 
condition h(m, 0) = 0, we get for m/0, 



h(m,t)= e -4- 2 ™ 2 r(t-0/L 2 f,{m,t')dt'. (A4) 
Jo 

Using Eq. (|A4|I and the properties of the noise, one can easily compute the two-point equal time correlation function 
and we get 



2D' L 2 
(h{m u t)h(m 2 ,t)) = — - 



TL 8w 2 mf 

This gives the two point correlation function in real space, 

2D' 



1 - e 



-8ir 2 mlrt/L 2 



(h(x 1 ,t)h(x 2 ,t)) = =^j- £ -±-, i-e-*^"/* 



%i+m2,0- 



(A5) 



(A6) 



Note that the sum in Eq. (|A6|I runs from m — — oo to m = oo but does not include the m = term, since we have 
used the fact that h{m = 0, t) — 0. In particular, putting x\ — x 2 — x, we get the on-site variance, which becomes 
independent of x as expected, due to the translational invariance. This gives 



<^ 2 (o)> = ^ E 11 



mi^O 1 



Y _ e -8ir 2 mlrt/L 2 



(A7) 
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In the stationary limit (t — > oo), one gets from Eq. (|A7|) 

<^(°)> = £ A" (A8) 

m— 1 

Using the identity, X)m=i = 1x2 l^i we S e t the formula for the on-site variance in the stationary limit 



{h2{0)) = & L ^ £ ^' (A9) 

where D = D' /2T. Since the Eq. is linear, it follows then that the single site height distribution is a pure Gaussian 
at all times. In particular, the stationary distribution is given by 

P st (ho) = 7= * . exp[-h 2 Q /2(h 2 (0))l (A10) 

where (h 2 (0)) is given in Eq. (|A9|) . 



